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Let {X{t) = {Xi{t), X 2 {t))'^, t £ R'^} be an R^-valued continuous locally stationary Gaussian 
random field with E[X(t)] = 0. For any compact sets Ai,A 2 C R^, precise asymptotic behavior 
of the excursion probability 

P ( max Xi (s) > u, max X 2 (t) > u], as u ^ 00 

y t^A.2 J 

is investigated by applying the double sum method. The explicit results depend not only on 
the smoothness parameters of the coordinate fields Xi and X 2 , but also on their maximum 
correlation p. 

Keywords: Bivariate Gaussian field, Bivariate Matern field, Excursion probability, Double ex¬ 
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1. Introduction 

For a real-valued Gaussian random field X = {X{t), t G T}, where T is the parameter 
set, defined on probability space {Ul, T, P), the excursion probability Pjsup^gy X (t) > u} 
has been studied extensively. Extending the seminal work of Pickands (1969), Piterbarg 
(1996) developed a systematic theory on asymptotics of the aforementioned excursion 
probability for a broad class of Gaussian random fields. Their method, which is called the 
double sum method, has been further extended by Chan and Lai (2006) to non-Gaussian 
random fields and, recently, by D^bicki et al. (2014) to a non-stationary Gaussian ran¬ 
dom field {X(s,t), (s,t) £ whose variance function attains its maximum on a finite 
number of disjoint line segments. For smooth Gaussian random fields, more accurate 
approximation results have been established by using integral and differential-geometric 
methods (see, e.g., Adler (2000), Adler and Taylor (2007), Azai's and Wschebor (2009) 
and the references therein). For Gaussian and asymptotically Gaussian random fields, 
the change of measure method was developed by Nardi, Siegmund and Yakir (2008) and 
Yakir (2013). Many of the results in the aforementioned references have found important 
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applications in statistics and other scientific areas. We refer to Adler, Taylor and Worsley 
(2012) and Yakir (2013) for further information. 

However, only a few authors have studied the excursion probability of multivariate 
random fields. Piterbarg and Stamatovich (2005) and D^bicki et al. (2010) established 
large deviation results for the excursion probability in multivariate case. Anshin (2006) 
obtained precise asymptotics for a special class of nonstationary bivariate Gaussian pro¬ 
cesses, under quite restrictive conditions. Hashorva and Ji (2014) recently derived precise 
asymptotics for the excursion probability of a bivariate fractional Brownian motion with 
constant cross correlation. The last two papers only consider multivariate processes on 
the real line R with specific cross dependence structures. Cheng and Xiao (2015) es¬ 
tablished a precise approximation to the excursion probability by using the mean Euler 
characteristics of the excursion set for a broad class of smooth bivariate Gaussian random 
fields on . In the present paper we investigate asymptotics of the excursion probability 
of non-smooth bivariate Gaussian random fields on R^, where the methods are totally 
different from the smooth case. 

Our work is also motivated by the recent increasing interest in using multivariate ran¬ 
dom fields for modeling multivariate measurements obtained at spatial locations (see, e.g., 
Gelfand et al. (2010), Wackernagel (2003)). Several classes of multivariate spatial models 
have been introduced by Gneiting, Kleiber and Schlather (2010), Apanasovich, Genton and Sun 
(2012) and Kleiber and Nychka (2012). We will show in Section 2 that the main results 
of this paper are applicable to bivariate Gaussian random fields with Matern cross¬ 
covariances introduced by Gneiting, Kleiber and Schlather (2010). Furthermore, we ex¬ 
pect that the excursion probabilities considered in this paper will have interesting sta¬ 
tistical applications. 

Let {X{t),t G R'^} be an R^-valued (not-necessarily stationary) Gaussian random 
field with E[A(t)] = 0. We write X{t) = {Xi{t), X 2 (t))'^ and define 

:=E[Xi(s)Xj{t)], i,j = l,2. (1.1) 

Let |t| := be the 1^-norm of a vector t G R^. Throughout this paper, we 

impose the following assumptions. 

i) rii(s, t) = 1 — Ci\t — s|“* -b o{\t — s|“’), where Ui G (0, 2) and > 0 (* = 1, 2) are 
constants. 

ii) \rii{s, t)| < 1 for all \t — s\ > 0, * = 1, 2. 

hi) ri 2 {s,t) = r 2 i{s,t) := r{\t — s|). Namely, the cross correlation is isotropic, 
iv) The function r(-) : [0, oo) —>• R attains maximum only at zero with r(0) = p G {0, 1), 
i.e., |r(t)| < p for all t > 0. Moreover, we assume ^'(O) = 0,r"(0) < 0 and there 
exists r] > 0, for any s G [ 0 , 77 ], r"(s) exists and continuous. 

The cross correlation defined here is meaningful and common in spatial statistics 
where it is usually assumed that the correlation decreases as the distance between two 
observations increases (see, e.g., Gelfand et al. (2010), Gneiting, Kleiber and Schlather 
(2010)). We only assume that the cross correlation is twice continuously differentiable 
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around the area where the maximum correlation is attained, which is a weaker assumption 
than that in Cheng and Xiao (2015) who considered smooth bivariate Gaussian fields. 

For any compact sets Ai,A 2 C we investigate the asymptotic behavior of the 
following excursion probability 

PI maxXi(s) > M, maxX 2 (t) > M ), as m —>■ oo. (1.2) 

y t^A.2 J 

The main results of this paper are Theorems 2.1 and 2.2 below, which demonstrate that 
the excursion probability (1.2) depends not only on the smoothness parameters of the 
coordinate fields Xi and X 2 , but also on their maximum correlation p. The proofs of 
our Theorems 2.1 and 2.2 will be based on the double sum method. Compared with 
the earlier works of Ladneva and Piterbarg (2000), Anshin (2006) and Hashorva and Ji 
(2014), the main difficulty in the present paper is that the correlation function of Xi 
and X 2 attains its maximum over the set D := {(s, s) : s £ Ai n A 2 } which may have 
different geometric configurations. Several non-trivial modifications for carrying out the 
arguments in the double sum method have to be made. 

This paper raises several open questions. First, the cases of ai = 2 or 02 = 2 have 
not been considered in this paper. The main difficulty is that, when (X\ = 2, the sample 
functions of Xi may either be differentiable or non-differentiable. In view of the method 
in this paper, the proof of Lemma 4.1 on the uniform convergence of finite dimensional 
distributions for bivariate process breaks down when ai = 2 or a 2 = 2. Studying the 
asymptotics of (1.2) when ai = 2 or/and a 2 = 2 requires different methods for dealing 
with differentiable or non-differentiable cases. When both Xi and X 2 have twice continu¬ 
ously differentiable sample functions, this problem has been studied by Cheng and Xiao 
(2015). The authors plan to study the remaining cases in their future work. Second, 
it would be interesting to study the excursion probabilities when {X{t), t £ M.^} is 
anisotropic or non-stationary, or taking values in with d > 3. In the last problem, the 
covariance and cross-covariance structures become more complicated. We expect that the 
pairwise maximum cross correlations and the size (e.g., the Lebesgue measure) of the set 
where all the pairwise cross correlations attain their maximum values (if not empty) will 
play an important role. 

The rest of the paper is organized as follows. Section 2 states the main theorems with 
some discussions and provides an application of the main theorems to the bivariate Gaus¬ 
sian fields with Matern cross-covariances introduced by Gneiting, Kleiber and Schlather 
(2010). We state the key lemmas and provide proofs of our main theorems in Section 3. 
The proofs of the lemmas are given in Section 4. 

We end the introduction with some notation. For any t £ , |t| denotes its ^^-norm. 

An integer vector k £ is written as k = (fci,..., fcjv). For k £ and T £ K.+ = [0, 00 ), 
we define the cube [kT, (k -|- 1)T] Ylf^iihT, {h + 1)T]. For any integer n, mesn(-) 
denotes the n-dimensional Lebesgue measure. An unspecified positive and finite constant 
will be denoted by Co. More specific constants are numbered by Ci, C 2 ,_ 
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2. Main Results and Discussions 


We recall the Pickands constant first (see, Pickands (1969); Piterbarg (1996); Dieker and Yakir 
(2014)). Let X = be a (rescaled) fractional Brownian motion with 

Hurst index a/2 G (0,1), which is a centered Gaussian field with covariance function 

]E[x(0x(s)] = i^r + kr -\t- si“- 

As in Ladneva and Piterbarg (2000) and Anshin (2006), we define for any compact 
sets §, T c , 

Ha{S,T) := [ e“-P(sup(x(s)-|s|“) >a;, sup(x(f)-|t|“) >a;)da;. (2.1) 

Jo ^ sGS tGT ^ 


Let Ha{T) = iLa(T,T). Then, the Pickands constant is defined as 


Ha ■= lim 

T—>-oo 


Ho.{[0,T]^) 

J'N 


( 2 . 2 ) 


which is positive and finite (cf. Piterbarg (1996)). 

Before moving to the tail probability of extremes of a bivariate Gaussian random field, 
let us consider the tail probability of a standard bivariate Gaussian vector (^, rf) with 
correlation p. It is known that (see, e.g., Ladneva and Piterbarg (2000)) 


P(^ > u,r] > u) =4>(m, p ){1 + o(l)), as M —>■ oo, 


where 


4'(u,p) := 


(1 + P)^ 


2ttu'^ 


exp - 


l + p 


The exponential part of the tail probability above is determined by the correlation p. 
As shown by Theorems 2.1 and 2.2 below, similar phenomenon also happens for the 
tail probability of double extremes of {X{t),t G K^}, where the exponential part is 
determined by the maximum cross correlation of the coordinate fields Xi and X 2 . 

We will study double extremes of X on the domain Ai x A 2 where Ai, A 2 are bounded 
Jordan measurable sets in . That is, the boundaries of Ai and A 2 have TV-dimensional 
Lebesgue measure 0 (see, e.g., Piterbarg (1996), p.l05). We only consider the case when 
Ai n A 2 ^ 0, in which the maximum cross correlation p can be attained. 

If mesN{Ai n A 2 ) ^ 0, we have the following theorem. 


Theorem 2.1. Let {X{t),t G K.'^} be a bivariate Gaussian random field that satisfies 
the assumptions in Section 1. If mesN{Ai C\ A 2 ) 0, then as u ^ 00 , 


PI maxAi(s) > u, maxA 2 (t) > u 

sGAi tGA 2 


= (27r) 2 (-r"(0)) 2 mesjv(Ai AA 2 )Ha^Ha^ 

X (1 -I- p){1 0 ( 1 )). 


(2.3) 
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If mesN{Ai C\ A 2 ) = 0, the above theorem is not informative. We have not been able 
to obtain a general explicit formula. Instead, we consider the special cases 

N N 

^1 = ^i,M X [Sj,Tj\ and A 2 = A 2 ^m x (2-4) 

j=M+l M+1 

where Ai^m and A 2 ^m are M dimensional Jordan sets with mesM{Ai^M H A 2 ,m) ^ 0 
and Sj < Tj < Rj,j = M + 1,...,N,0< M < N — 1. For simplicity of notation, let 
meso(-) = 1. Our next theorem shows that the excursion probability is smaller than that 
in (2.3) by a factor of . 

Theorem 2.2. Let {X{t),t £ R.'^} be a bivariate Gaussian random field that satisfies 
the assumptions in Section 1, and let Ai,A 2 be as in (2.4) sueh that mesM{Ai^M O 
A 2 ,m) > 0. Then as u ^ 00 , 


IX, liid^ . 

t£A2 


PI maxJfi(s) 

\ s^A\ 

N N 

= (27r)^(-r"(0))-^cr 


X 


(1+p) 


2N 


II '.'1 ^2 (^l,M O j42,M ) 

M+N{^ + ^ — 2),j,f w. , ..XX 

“1 “2 U ^“1 “2 p)(l + o(l)). 


(2.5) 


Remark 2.1. The following are some additional remarks about Theorems 2.1 and 2.2. 

• The excursion probability in (1.2) depends on the region where the maximum cross 
correlation is attained. In our setting, the maximum cross correlation p is attained 
on D := {(s, s) | s £ Ai n 242 }. 

• For Theorem 2.2, let us consider the extreme case when M = 0, i.e., Ai fl .42 = 

{(Ti,..., Tjv)}. The exponential part still reaches although the maximum 

cross correlation p is attained at a single point. 

• To compare our results with Anshin (2006), we consider a centered Gaussian process 

{X{t) = {Xi{t),X 2 {t)),t £ R} and Ai = A 2 = [0,T]. In our setting, the cross 
correlation attains its maximum on the line D = {(s, s) | s £ [0,T]}, while in 
Anshin (2006) it only attains at a unique point in [0,T] x [0,T] because of the 
assumption C2. This is the reason why the power of u in our settings is ^ ^ ~ 3 

instead of + ^ — 4 in Anshin (2006). 

• Even though Theorem 2.2 only deals with a special case of Ai, A 2 with mesN{Ai fl 
A 2 ) = 0, its method of proof can be applied to more general cases provided some 
information on Ai and A 2 is specified. The key step is to reevaluate the infinite 
series in Lemma 3.5. 


We end this section with an application of Theorems 2.1 and 2.2 to bivariate Gaussian 
random fields with the Matern correlation functions introduced by Gneiting, Kleiber and Schlather 
( 2010 ). 
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The Matern correlation function M{h\v, a), where a > 0,1^ > 0 are scale and smooth¬ 
ness parameters, is widely used to model covariance structures in spatial statistics. It is 
defined as ^ 

M{h\v,a) := ^^{a\h\yK^{a\h\), (2.6) 


where Ki, is a modified Bessel function of the second kind. In Gneiting, Kleiber and Schlather 
(2010), the authors introduce the full bivariate Matern field df(s) = (Xi(s),^ 2 ( 5 )), i.e., 
an M^-valued Gaussian random field on with zero mean and matrix-valued covariance 
functions: 


where Cij{h) := 


C{h)-- 

( C^^[h) Ci2{h) \ 

V C2l{h) C22{h) ) ’ 

(2.7) 

]E[Ai(s + /i)Aj(s)] 

are specified by 


Cii{h) = 

alM{h\vi,ai), 

(2.8) 

C22{h) = 

alM{h\v2,a2), 

(2.9) 

Ci2{h) = 

C2i{h) = paia2M{h\vi2,ai2). 

(2.10) 


According to Gneiting, Kleiber and Schlather (2010), the above model is valid if and 
only if 


2 / r(^i + N/2)Tiu2 + N/2) V{u^2? 

^ - n^i)V{y2) r(^i2 + IV/2)2 

f^2 I ^2\2izi2+IV 

X inf_ *-”12 + ^ ) _ 

t >0 (af + t 2 )>'l-HJV/ 2 (a 2 i 2 ^^ 2 +IV /2 ■ 

Especially, when oi = 02 = 012 , condition (2.11) is reduced to 

2 ^ r(i2i + N/2)r(iy2 + N/2) r{i^i 2 )^ 

^ - r(i.i)r(zz2) r(j.i2 + fV/2)2’ 


( 2 . 11 ) 


( 2 . 12 ) 


in which case the choice of p is fairly flexible. 

Here we focus on a standardized bivariate Matern field, that is, we assume cti = 172 = Ij 
Or = 02 = 012 = 1 and p > 0. Moreover, we assume 1 ^ 1 , V 2 G (0,1) and 1^12 > 1. In this 
case, the bivariate Matern field {X(t), t G R^} satisfies the assumptions in Section 1. 
Indeed, Assumption i) in Section 1 is satisfied since 

M{h\i'i,a) = 1 - + odtp'"'). 


where ci = 2 :i^ir(i+!-) ’ * = 1)2 (see, e.g., Stein (1999), p. 32). Assumption ii) holds imme¬ 
diately if we use the following integral representation of M {h\v^ a) (see, e.g., Abramowitz and Stegun 
(1972), Section 9.6) 


M{h\h', a) 


2T{i^+l/2) 

0Fr(i/) 


cos(a|h|r) , 

_ ffr 

(^l^j.2y+i/2 “ ■ 


0 


(2.13) 
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Assumption iii) holds by the definition of cross correlation in (2.10). For Assumption iv), 
we only need to check the smoothness of M{h\i^,a). By another integral representation 
of M{h\i',a) (see, e.g., Abramowitz and Stegun (1972), Section 9.6), i.e., 

= Tr+i/2 )r() r 


one can verify that M{h\v, a) is infinitely differentiable when \h\ ^ 0. Meanwhile, M"(0|zz, a) 
exists and is continuous when v > 1 which can be proven by taking twice derivatives to 
the integral representation in (2.13) w.r.t. \h\. So Assumption iv) holds. 

Applying Theorem 2.1 to the double excursion probability of Ar(s) over [0,1]^, we 
have 


( max Xi(s) > u, max X 2 {t) > u] 
VsGfo.ilJv tefo.il" / 


P 

sG[0.1]JV ' ' ' tG[0.1]^ 

= (27r)^(-C" (0))-^Ci^C2^(l+ + 

p)(i-p o(l)), as u —>■ oo. 


Secondly, when the two measurements are observed on two regions which only share 
part of boundaries, we use Theorem 2.2 to obtain the excursion probability. For example, 
if A'i(s) are observed on the region [0,1]^ and ^ 2 ( 5 ) on [0, x [1, 2], then as u —>■ 00 , 


max A'i(s) > u, max A' 2 (t) > u 

iG[0,l]'^ tG[0,l]"‘-ix[1.2] 


= ( 2 ^)^(-cr 2 (o))- 


T 21^2 


(1 + p) 




H2i'i H2V2 


Xy^(«'i + ^2 ^'I'(u, p)(l + 0(1)). 


3. Proofs of the main results 

The proofs of Theorems 2.1 and 2.2 are based on the double sum method (Piterbarg 
(1996)) and the work of Ladneva and Piterbarg (2000). Since the latter deals with the 
tail probability P(maxtg[j’j X{t) > u, maxtgjj’j T'^] X{t) > u) oi a univariate Gaussian 
process {X{t),t S R}, their method is not sufficient for carrying out the double sum 
method for a bivariate random field. 

Lemmas 3.1 and 3.2 below extend Lemma 1 and Lemma 9 in Ladneva and Piterbarg 
(2000) to the bivariate random field {(Xi(t),X 2 (t)), t G R^}. Moreover, we have strength¬ 
ened the conclusions by showing that the convergence is uniform in certain sense. This 
will be useful for dealing with sums of local approximations around the regions where 
the maximum cross correlation is attained. The details will be illustrated in the proof 
of Theorem 2.1 (see, e.g., (3.10), (3.21)). In the following lemmas, {X{t)^ t G R'^} is a 
bivariate Gaussian random field as defined in Section 1. 
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Lemma 3.1. Let Su and be two -valued functions of u and let t„ '■= tu — Su- For 
any compact sets S and T in , we have 


PI max -^i(s) > max X 2 {t) > u 

^ s£Su-\-u~^^°‘iS 


ji + pf 


--He 


X u ^exp 


.(1 + P)^ 
..2 


Hr. 




,(1 + P)^ 


(3.1) 


(1 + 0 ( 1 )), 


1 + ^(kul), 

where o(l) —>■ 0 uniformly w.r.t. t„ satisfying |tk| < Coy/logu/u as u —>■ oo. 

Lemma 3.2. Let s^, tu and Tu be the same as in Lemma 3.1. For all T > 0, m, n G , 

we have 


max -^i(s) > w, max ^ 2 (t) > u, 

ses„+u-2/“i [O.T]'^ tGtu+«“^/“2[0,T]^ 

max -^ 1 ( 5 ) > u, max ^ 2 (t) > u 

[mT,(m+l)T] [nT,(n+l)T] 


(i + p)^ 


27ry^l — y2 


e T+HFUTT He 


cY°‘^[0,T]^ [mT, (m + 1)T] 


(1 + P) 

^1/02 rp, rjTiN „l/a 2 


(1 + p)^ 


X He 


(3.2, 

V (1 + p )“2 (1 + p )“2 / 


where Ha{-,-) is defined in (2.1) and o(l) —>■ 0 uniformly for all s„ and tu that satisfy 
\tu\ + CoVlogw/w as u ^ 00 . 

Proofs of Lemmas 3.1 and 3.2 will be given in Section 4. Now we proceed to prove our 
main theorems. 

Proof of Theorem 2.1. Let 11 = Ai x >12, 5{u) = C-v/log u/u, where C is a constant 
whose value will be determined later. Let 


V = {(s,t) GlI: |t — s| < S{u)}. 


(3.3) 


Since 


I J {-^i(s) > u,X2{t) > u)} ) < P( maxXi(s) > u, ma.xX2{t) > u] 

'' ‘eA2 J 

+ 1 ?^ U {-’^i(s) > '“,-’^2(0 > m)}^ + IJ {Xi(s) > M,X2(t) > n)}^, 


' {s,t)ev 


' {s,t)en\v 

















Tail Asymptotics for Extremes of Bivariate Gaussian Fields 


9 


it is sufficient to prove that, by choosing appropriate constant C, we have 


and 


[J {Xi{s) > u, X 2 it) > u)} 

(s,i) 

= ( 27 r)'^(—(1 + p)~^^'^~^'^~^'^mesN{Ai n A2) 

X Ha^Ha^ p)(l + 0(1)), as u —>■ 00 


IP (U(s,t)en\i>{^i('S) > u,X2{f) > u)}) 
' (U(s,t)Gi’{^i(s) > u,X2{t) > w)}) 


lim 

U—^OO 


= 0 . 


(3.4) 


(3.5) 


We prove (3.4) first. For any fixed T > 0 and * = 1,2, let di{u) = Tu and, for any 
k = (fci,..., fc^r) G define 


N 


{kj + l)di{u)] = [kdi{u), (k + l)di{u)]. 
i=i 


(3.6) 


Let 

C = {(k, 1) : X n P 7 ^ 0} and C° = {(k, 1) : A^ x Ap^ C V}. (3.7) 

It is easy to see that 

U AP’xAp^CPC y Ap^xAp). 

(k,I)GC° (k,l)GC 


Thus the LHS of (3.4) is bounded above by 


Let 



U,X2{t) > It)} 


< E p 

(k,l)GC 


max Ai(s) > u, max X 2 (t) > u 
sGA^i> tGAp) 


= E p 

(k,l)GC 


max Aii(s) > it, max X 2 {f) > u 

sGkcii(u)+A(|'‘^ tGld2(ti)+A(^^^ 


Tki := ld2{u) — k(ii(it) 

= ihd2{u) - kidiiu), ..., lNd2{u) - k^di^u)). 


(3.8) 


(3.9) 
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For (k, 1) e C, |Tki| < S{u) + VN{di (m) + d 2 {u)) < 26{u) for all u large enough, since 
di(u) = o(6(u)) and d 2 (u) = o(6(u)), as u —> oo. By applying Lemma 3.1 to the RHS of 
(3.8), we obtain 



,X 2 {t) > u)} 


^ (l + p)^(l + y(n)) ^ ^ f 

- 2ny^u^ “^1, (1+p)* ) “=1, (1 + p)* ) 

X V expf- ^ \ 

m ^ 1 + Kkkil); 


= H„ 


1 1 2 

V ( 1 +P)- . 

X ^ exp -i 
(k.l)GC ^ 


Hr. 


, (1+p)^ , 

1 


^'(■u,p)(l + 7(u)) 


l+J'dTkll) 1+P 


(3.10) 


where the global error function "f{u) —>■ 0, as u —>■ c». The uniform convergence of (3.1) in 
Lemma 3.1 guarantees that the local error term o(l) for each pair (k, 1) G C is uniformly 
bounded by 'y{u). 

The series in the last equality of (3.10) is dealt by the following key lemma, which 
gives the power of the threshold u in (3.4). 


Lemma 3.3. Recall the set C defined in (3.7). Let 

Then, under the assumptions of Theorem 2.1, we have 

h{u) = (27r)'^/2(-r"(0))-'^/2(l + p)^mesN{Ai n A 2 ) 

X y^(~+sy-i)(i _|_ q(i)), fj:g y (3.12) 

Moreover, if we replace C in (3.11) by C° defined in (3.7), then (3.12) still holds. 


We defer the proof of Lemma 3.3 to Section 4 and continue with the proof of Theorem 
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2.1. Applying (3.12) to (3.10), we obtain 
P(^ IJ {Xi{s)>u,X 2 it)>u)} 


' {s,t)ev 


< {27r)^{-r''iO))-^il + py^T-^‘^mesNiA,nA2)H^, 


, (1+p)^ , 


X 




+ p)(l + 7 i(u)), 


V (l + p)“2 
where 71 (u) —>■ 0, as u —>■ 00 . Hence, 

' (U(s,t)GX>i^i(s) > u,X2{t) > m)}) 


lim sup ■ 

U—XX) 




< ( 27 r) 2 (-r"( 0 )) 2 (1 + p)^mesNiAi n A2) 

..r-2Nrr /cy“^[0,r]^' 

X J. irLf^-y \ ry I JEL 0^2 \ O 


(1+p)^ 


(1+p)^ 


(3.13) 


(3.14) 


The above inequality holds for every T > 0. Therefore, letting T —>■ oo, we have 
' (U(s.t)Gi’{^i(s) > u,X2{t) > u)}) 


lim sup ■ 

U—¥(X) 


< (2^)^(-r"(0))-^ 


X C2 ^ (1 + p)~ ^'^~'~^~^’mesN{Ai n A2)Ha-^Ha2- (3.15) 

On the other hand, the lower bound for LHS of (3.4) can be derived as follows. Let 

H = {(k, 1, k', 1') : (k, 1) ^ (k', 1'), (k, 1), (k', 1') G C}. (3.16) 


By Bonferroni’s inequality and symmetric property of B, the LHS of (3.4) is bounded 
below by 



,X2it) > U} 


> E p 

(k.l)GC° 


max Xi(s) > u, max X 2 (f) > u 

seAy> teAy> 


1 

2 


E 

(k,l.k',l')GB 


max Xi(s) > u, max X 2 (t) > u, 

tGA® 


max Ai(s) > u, max X 2 {t) > u 
seA<^V ieA<?) 


(3.17) 
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= Si - S 2 . 


Since C° and C are almost the same, a similar argument as in (3.10)^ (3.15) shows that 
El is bounded below by 


El > {27T)-{-r"iO))--{l + py^mesNiAinA2)T-^^^Ha,, 


-2N^ 


^N 


X i/o 




(1 + p)‘ 


V (i + p )“2 / 

where 72 (m) —?> 0, as m —>■ 00 . Hence, letting T —>■ 00 , we have 

El 


lim inf ■ 




>( 2 ^)^(-r"( 0 ))-^Ci-cr 


x(l + p) ^'>mesN{Air\ A 2 )Ha^Ha^. 

Next, we consider S 2 in (3.17). To simplify the notation, we let 


(3.18) 


(3.19) 


/(k, 1 , k', 1 ') := P( max Xi(s) > M, max X 2 (t) > u, 

tGAp> 


max Xi(s) > It, max X 2 {t) > u 
sGA^V * 6 A(f’ 

k' 1' 


For m = (mi,..., ttin) S , let 


Ha,e(m) 4 He 


c^/“[ 0 ,T]^ ci/“[mr, (m+l)r] 


(1 + p)- (i + p)- 

Rewriting E 2 and applying Lemma 3.2, we obtain 

^’ = )EfE + E + E )t(i‘,i,kM' 


(3.20) 


(k,l)GC ^ (k'.l')GC (k',l')GC (k',l')GC 

k'=kaVl kVk,l'=l kVklVl 

.(l+p)^(l+73(«)) ^ 


47r\/l-p2ii2 


l + r(|-rj^ll) 


-He^cAO) ^- 2 .C 2 ( 1 '- 1 ) 


(k',l')GC 

k'=k,lVl 


+ ne,,cA0) ^ai.ci(k'-k)+ ^ H„,,,,(k'-k)H„,.c,(l'-l) 


(k',l')GC 

kVk.l'=l 


< 


(1 +p)^(l + 73 ( 11 )) 

Yn\J\ — 


^ ) g rnTFuaTT 
(k,l)GC 


(k'.l')GC 

k+k.l+l 

^ai,ci( 0 ) ^ "Hc^.csCn) 

ni^O 
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“t" ^a2jC2(0) ^ ^ ^ ^ ^ai ,ci (^^)^CK2 ,C2 (^) j ; 

m/0 m/0,n^0 


(3.21) 


where 73 ( 14 ) —> 0, as it —>■ cxd. According to the uniform convergence of (3.2), the local 
error term o(l) for each pair (k', 1') G C is bounded above by 73(14) . To estimate 'Ha,c{-), 
we make use of the following lemma, whose proof is again postponed to Section 4. 

Lemma 3.4. Recall 'Ha,c{') defined in (3.20). Let iq = argiiiaX]^<i<Ar|iiii|- Then there 
exist positive constants Ci and Tq such that for all T > Tq, 


Consequently, 


■Hc.A0)<CiT^; 

'Ha,ci^) < CiT^~^, when \mA = 1 ; 

■Ha,c(m) < when |rn,J > 2. 


^ H„.c(m) < 

mez"\{0} 


(3.22) 

(3.23) 

(3.24) 

(3.25) 


Applying Lemmas 3.3 and 3.4 to the RHS of (3.21), we obtain 
S2 < C'0(l + Pg^ 73(14)) ^ 24 V -1 ^ 


1 + l’(l'Tcll) 


4’r\/l-p'i4" (k,i) 6 C 

< Co( 27 r)'^(-r"( 0 ))“'^(l + p)^mesjv(Ai n A 2 )T~^ 

X p){l + 74 ( 14 )), 

where 74 ( 14 ) —>■ 0, as it —>■ 00 . By letting it —>■ 00 and T —>■ 00 successively, we have 

S 2 


lim sup ■ 




r-i). 


= 0 . 


U-S-OO ' 02 "'’^(u,p) 

By combining (3.17), (3.19) and (3.27), we have 


(3.26) 

(3.27) 


lim inf ■ 

U—^OO 


(U(s,t)Gl5{-^l(s) > U,X2{t) > U)}) 


> lim inf ■ 


Si 


4'(u,p) 

— lim sup ■ 


S 2 


^^4>(i4,p) U^OO y^(oi+o2 ^^4>(lt, p) 

> (27r)^(-r"(0))“'^Ci“^C2“^ (1 + p)“^^^'^^“^^TOesiv(Ai n A2)77aii7Q 

It is now clear that (3.4) follows from (3.15) and (3.28). 

Now we prove (3.5). Define 


(3.28) 


Y{s,t) := Ai(s) +X 2 (t), for (s,t) G 11 \ 27. 


(3.29) 
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For X = {si,ti),y = ( 52 ,^ 2 ) G H \ X>, let \x - y\ = Then we can 

verify that 

E|F(cr)-y(y)p<Co|x-yr“(“^’“^\ yx,y€U\V. (3.30) 

By applying Theorem 8.1 in Piterbarg (1996), we obtain that the numerator of (3.5) is 
bounded above by 



< P 


f max Y(s, t) > 2u 
V(s,t)Gn\-D 


< Cou 


exp 


1 + max(s^t)gn\D r{\t - s\) J ' 


(3.31) 


Since r{\t — s|) = p + ir"(0)|t — sp(l + o(l)) and r(-) attains maximum only at zero, we 
have 


max r{\t — s|) < p 


^{-r"{0))d^{u) 


(3.32) 


for u large enough. So (3.31) is at most 


I I 2JV 

CqU min{cji,02) eXp I — 


1 + p - |(-r"(0))(52(u) 


I I 21V 

< CqU ^min(ai,Q2) exp ( — 


1 +P 


exp 


O + P? 


27r-\/l - p^Co 1 + 
(1 + p)^ "" 


-r"(0) ^2 

'I'(m,p), 


(3.33) 


where the inequality holds since ^ + ^,Vx > y. Compare (3.33) with (3.4), it is 

easy to see (3.5) holds if and only if 


2N -r"(0) 

1 H- —— 

min(Q;i,a2) 3(1 + p)^ 

Hence, by choosing the constant C satisfying 

2 

min(ai, a 2 ) 


C > 


3(1 + p)^ 

-r"(0) 


("( 


<N 






(3.34) 


(3.35) 


we conclude (3.5). 


□ 


Proof of Theorem 2.2. From the proof of Theorem 2.1, we see that the exponential 
decaying rate of the excursion probability is only determined by the region where the 
maximum cross correlation is attained. In the case of mesN{Air\A 2 ) = 0 but Hi(" 1^2 ^ 0, 

the exponential part, e~i+p, remains the same. Yet, the dimension reduction of Hi n H 2 
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does affect the polynomial power of the excursion probability, which is determined by 
the quantity 

in Lemma 3.3. Under the assumptions of Theorem 2.2, the set C and the behavior of h{u) 
change. We will make use of the following lemma which plays the role of Lemma 3.3. 

Lemma 3.5. Under the assumptions of Theorem 2.2, we have 

h{u) = (27r)^/^(-r"(0))“/^"^(l + p)^^"“T"2^TOesM(^i,M n A 2 ,m) 

X ^M+Ar(—+—- 2 ) (1 ^( 1 )), y (3.36) 

Moreover, if we replace C with C° defined in (3.7), then the above statement still holds. 

The rest of the proof of Theorem 2.2 is the same as that of Theorem 2.1 and it is 
omitted here. □ 


4. Proof of Lemmas 


For proving Lemma 3.1, we will make use of the following 


Lemma 4.1. Let Su and he two -valued functions of u and let t„ := tu — Su- For 
any compact rectangles S and T in define 

fu{s) := u(Xi(su + - u) X, Vs e S, 

r]u{t) := u(X2(tu + - u) + y, Vt G T (4.1) 

and for any t G let 

m ■■= V^iXi(t) - v(t) ■■= V^X2{t) - (4.2) 

where Xi(^) 1 X 2(0 independent fractional Brownian motions with indices ai/2 

aud 012 / 2 , respectively. Then, the finite dimensional distributions (abbr. f.d.d.) of {^ui'), Vui')), 
given Xi(su) = tt — X 2 (tu) = w — converge uniformly to the f.d.d. of (^(•), r](-)) for 
all Su and tu that satisfy Ir^l < Co\t\ogu/u. Furthermore, as u ^ 00 , 


X y 

max^„(s) > x,ma,xriu(t) > y Xi(s„) =u -,X 2 (<„) =u - 

sGS iGT U U 


maxf(s) > x,maxu(t) > y 
sgS tel 


(4.3) 


where the convergence is uniform for all Su and tu that satisfy \tu\ < Cov/logu/u. 
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Proof. First, we prove the uniform convergence of finite dimensional distributions. Given 
=u—^, X^itu) = the distribution of the bivariate random field 
is still Gaussian. Thanks to the following lemma (whose proof will be given at the end 
of this section), it suffices to prove the uniform convergence of conditional mean and 
conditional variance. 


Lemma 4.2. Let X{u,Tu) = t„), ..., X„(m, t„))^ be a Gaussian random vec¬ 
tor with mean /i(M,T„) = Tu), ...,/x„(u,Tii)^ and covariance matrix Y{u,Tu) with 

entries aij{u,Tu) = Cov{Xi{u, Tu), Xj{u,Tu)), i,j = 1,2,...,n. Similarly, let X = 
(Xi,...,X„)^ be a Gaussian random vector with mean /i = (/xi,...,/i„) and covari¬ 
ance matrix E = Assume that E is non-singular. Let F'«(-) and F{-) be the 

distribution functions of X(u,Tu) and X respectively. If 

lim max |/x, (u, Tu) —/x, | = 0, 

U—>-(X> Tu 

lim max|tTi, (u,r„) - tTi, | = 0, *, j = 1, 2,..., n, (4.4) 

U—^OO Tu 

then for any x G , 

lim max|Fu(a::) — F{x)\ = 0. (4.5) 

U—¥00 Tu 


We continue with the proof of Lemma 4.1 and postpone the proof of Lemma 4.2 to 
the end of this section. Recall that, for two random vectors X,Y G K™, their covariance 
is defined as Gov(X,E) := E[(X — EX)(E — EE)^] and the variance matrix of X is 
defined as Var(X) := Gov(X, X). The conditional mean of (^u(t), r]u{t))^ given Xi(su) = 
■u - f, X 2 itu) = M - is 


E 


Ut) 

Vuit) 

+ Gov 


Xi(s„) = u - f 

X2itu) =U-\ 


= E 


iu{t) 

Vu{t) 


XliSu) 

X2(tu) 


Ut) 

Vu{t) 


Var 


Xiisu) 

X2{tu) 


— U + X 

—u^ + y 

1 

-r(|r, 


+ 


_ rii{su+u s„) 

V r{\Tu + 

u-^ 


-npu 

1 


U — 


ai{u) 
02 (m) 



r{\Tu-u \ 

r22(tu + u-'^/°^H,tu) j 


(4.6) 


where 


ai{u) 


u'^{l-rii{su + u '^/°‘H,Su))-u'^{r{\Tu-u - r(|r„|)) 

1 + ^(k«l) 


{x - yr{\Tu\)){l - rii{su + u ‘^^°H,Su)) 
l-r‘^{\Tu\) 
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(y - xr{\Tu\)){r{\Tu\) - r{\Tu - u 


(4.7) 


and 


, . U^il-r 22 {tu + u - u^(r(|r„ + u - r(|T„|)) 

=-rT7(ra)- 

{y - xr{\Tu\)){l - r22{tu + u~'^/°^^t,tu)) 
l-r2(|r„|) 

jx - yri\Tu\))ir{\T^\) - r{\Tu + 

l-rH\ru\) 

Applying the mean value theorem twice, we see that for u large enough, 


(4.8) 


\r{\Tu+u ^^“t|) - r(|ru|)| < |u • max 

s is between 
|ru| and 

< • ma x k^(s)| 

|s|<2Co y/Togu/u 


< \u • max_ ( |s| • max |r"(t)| 

\s\<2CoVt^/u \ |t|<|d 


r'(s)| 


< 2C'o|t|\/logu • u ^ • m ax W{t)\ 

|t|<2Co y/logu/u 

< 4C'o|r"(0)||<|\/log'u • 14“^“^/“, 


(4.9) 
oo and 


where the second inequality holds because of = o(Vlog u/u), as u 

the last inequality holds since r"(-) is continuous in a neighborhood of zero. Thus (4.9) 
implies that, as u —>■ oo. 


u'^\r{\Tu+u ^/“t|) - r(|ru|)| < 4C'o|r"(0)||t| Vlogu • 0, 


(4.10) 


where the convergence is uniform for all s„ and t^ that satisfy |tii| < CoVlog u/u. We 
also notice that for i = 1,2 and all s G K.^, 


l — rii(s + u ^^°‘t,s)=CiU ^|t|“’+o(u ^),asu—>- 00 . 
By (4.6), (4.10), and (4.11), we conclude that, as u —>■ oo. 


E 


Ut) 

Vu{t) 


Xl{Su) =U - ^ 
X2{ty,) = U - I 


ci|*l° 


1 +P 

C 2 itr" 

i+p 


(4.11) 


(4.12) 


where the convergence is uniform w.r.t. Su and tu satisfying |r„| < Coy/logu/u. 

Next, we consider the conditional covariance matrix of (fuit) — fu(s), yu(t) — 77„(s))^. 


Var ( ( 


Ai(s„) \ 

X2(tu) j 
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= Var 


X Var 


Vuit) - 7 ?„(s) 

X2{tu) 


— Cov 


Cov 


Vuit) - r]u{s) 

Ut)-Us) 

riu{t) - r]u{s) 


XliSu) 

X2{tu) 

XliSu) 

X2itu) 


(4.13) 


Let hu{t,s) := r{\Tu + u — u Applying (4.10) and (4.11), we obtain 


Var 


Ut)-Us) \ 

Vu{t)-Vu{s) ) 


( 2u^{l-rii{su+ 

y-2/aig^ + rt“^/“it)) 

U^{hu{t,t) - hu{s,t) 

V 4 “ '^)) 


u^{hu{t,t) - hu{s,t) \ 

4“ ^ui.^1 '^)) 

2u^{\ - r22{tu+ 

tu + M“^/“^<)) / 


/ 2 ci|t-s|“ni + o(l)) o(l) 

V o(l) 2c2|t-s|“=(l + o(l)) 


(4.14) 


where o(l) converges to zero uniformly w.r.t. Tu satisfying |ru| < C^y/log u/u, as m —>■ oo. 
Also, we have 


( Ut)-Us) 

V Vu{t) - Vu{s) 


Ai(sJ \ 

X2itu) j 


/ u{rii{su + u s„) u(r(|r„-u \ 

-rii{su + u“^/“is,s„)) -r{\Tu - m“^/“is|)) 

u(r(|r„ + u“^/“=t|) u{r22{tu + <„) 

\ -r(|r„ + u“ 2 /“=s|)) -r22{tu + u~‘^^°“^s,tu)) J 

_ ( 0 ( 1 ) 0 ( 1 ) \ 

- 1 , 0 ( 1 ) 0 ( 1 ) ) ’ 

as u —>■ 00, and 


Var 


Ai(s„) 

X2{tu) 



1 

l-r 2 (|r„|) 


1 -r{K\) 

-r{\Tu\) 1 


By (4.13) - (4.16), we conclude that as u —>■ oo. 


f f ^u{t) - ?„(s) 
V V Vu{t) - vu{s) 


XliSu) 

X2{tu) 


2 ci\t — s|“i 
0 


0 ^ 
2 C2 |t-s|“= J ’ 


(4.15) 


(4.16) 


(4.17) 


where the convergence is uniform w.r.t. satisfying \tu\ < Cov/log u/u. Hence, the 
uniform convergence of f.d.d. in Lemma 4.1 follows from (4.12), (4.17) and Lemma 4.2. 

Now we prove the second part of Lemma 4.1. The continuous mapping theorem (see, 
e.g., Billingsley (1968), p. 30) can be used to prove (4.3) holds when and tu are fixed. 












Tail Asymptotics for Extremes of Bivariate Gaussian Fields 


19 


Since we need to prove uniform convergence w.r.t. Su and tu, we use a discretization 
method instead. Let 


f{u,x,y) 


:= PI max£„(s) > x, maxr7„(t) > y 

' sGS tGT 


=u - -,X 2 {tu) =u-- 
u u 


and 


f{x,y) 


maxf(s) > X, maxTiit) > y 
sgS teT ^ ' 


(4.18) 


(4.19) 


Without loss of generality, suppose that S = [a, b]^ and T = [c, where a < b,c < d. 
For any <5 G (0,1), let m = , n = and let 


Sra ■= {sk I Sk = {Xki, Xkt^), k = (fci, ..., ) G {0, 1, ...,771+ 1}^}, 

Tn ■= {ti I h = {yi^, ..., yi^), 1= {Ii,...,In) G {0,1,...,n + 1}^}, 

where Xi , yi are defined as 

a = Xq < Xi < ■ ■ ■ < Xm < Xm+l = b, Xi = a + iS, 7 = 0, 1, . . . , 777, 
c = 7/0 < 2/1 < • • • < yn < 7/n+l = d, ?/i = C + id, 7 = 0, 1, . . . , 77. (4.20) 

Then [a, 6]^ x [c, d]^ can be divided into d-cubes with vertices in Sm x Tn- 
The function f{u,x,y) in (4.18) is bounded below by 

fm,n{u,x,y) := FI max (s) > x, max pu (t) > y 

Y S^Sm i^Tn 

Xi{su)=u--,X2{tu) = u-'^) (4.21) 

u u J 

and is bounded above by gm,n{u, x, y) which is defined as 


max fu{s) > X — e, ma.xr]u{t) > y — e Xi(su) =u - , W 2 (t„) = u - 

seSm teTn u u 

( X y 

max ^ii(s) > a;, max 5«(s) < a; — e Xi(su) =77-,X 2 (t„) = u - 

sGS sGSm 77 77 

( X y 

max77„(t) > y, max7y„(t) <y - e Xi(s„) =u -,^2(tii) =77 - 

t^X'n U U 

= fm,niu, X - e, 2/ - e) + Sm,niu, X, y) + tm,n{u, X, ?/), 


(4.22) 


where e > 0 is any small constant. Let 


fm.nix^ 7 /) 


maxf(s) > X, max 77 ( 2 ) > y 

sG 5„ tGTn 


(4.23) 
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Since the finite dimensional distributions of ^«(‘)) converge uniformly to those of 

(?(•)>»?(•))> we have 


lim max \fm,n{'^,x,y) - fm,n{^,y)\ = 0- 

u—foo |T-u|<Co>/logli/'U 

The continuity of the trajectory of yields 

n—>-oo 

By (4.24) and (4.25), we conclude 

lim lim max |/,„,„(u, a;, y) -/(x, y)| = 0. 
lT„|<co^/I^/■u 

Let us consider the conditional probability Sm,n{u,x,y) in (4.22). 
y') 

<P( max \^u{s) - ^uit)\ > e Xi{su) = u - X 2 itu) = u -- 
y|s-i|<(5 u u 

< ( max IC^i(s) - Cu(i)| Xi{su) = u - X 2 {tu) = u -- 

e vb“*l <'5 u u 

1 


(4.24) 


(4.25) 


(4.26) 


=-Ep„ I max |x(s) — x(t)| ), 

|s — 1|<(5 


(4.27) 


where is the probability measure on ((7(8),^((^(S)) defined as 

P„(A) := P(^C«(-) G ^ 7fi(sti) =u - -,X2{tu) =u-- 
\ u u 

for all A € B(C(S)) and x(-) is the coordinate random element on ((7(§),S((7(S)), P„), 
i.e., x{t,uj) = u}{t), 'ioj G (7(8) and t G 8. Consider the canonical metric 

d„(s, t) : = [Ep„ (|x(s) - x(t)p)] 

1^2 

= E(|e.(s)-eu(0l" Xi{Su)=U--,X2{tu) = U-^)\ . 

L V u U / J 

By (4.17), for u large enough and all ««, tu such that \tu\ < (7o-\/logu/u, we have 

du(s,t) < 2ycr|s - (4.28) 

which implies Vs G 8, 

{t G 8 I |t - s| < (e/2ycr)~ } C {t G 8 I d„(s,t) < e}. 

Hence 


iVd„(S,e) 


(4.29) 
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where (S, e) denotes the minimum number of d^-balls with radius e that are needed 
to cover S. By Dudley’s Theorem (see, e.g., Theorem 1.3.3 in Adler and Taylor (2007)) 
and (4.28), we have 

/ \ _ 

Ep^ ( max |a;(s) — a;(t)| ) <K ^/\og Nd^ (S, e) de, (4.30) 

/ Jo 

where AT < oo is a constant (which does not depend on 6) and, thanks to (4.29), the last 
integral goes to 0 as 5 —>■ 0 (or, equivalently, as to —>■ oo, n —>■ oo). By (4.27) and (4.30), 
we conclude that 


liin^limsup m ax \sm,n{u,x,y)\ = 0. 

•OO \Tu I u\/u 


m—>-oo 
n—>-oo 


(4.31) 


A similar argument shows that 


lim^limsup m ax \tm,n{u,x,y)\ = 0. 

oo \Tu\<Coy/\og u\ / U 


m—>-oo 
n—^oo 


(4.32) 


Since 


\f{u,x,y) - f{x,y)\ < \fmA'>^,x,y) - f{x,y)\ + |5m,„(M,a;,y) - /(a;,y)| 

< |/m.n(w,a;,y) - /(a;,y)| + \ fm,n{u,x- e,y - e) - f{x - e,y - e)\ 

+ \ f{x - e,y - e) - f{x,y)\ + \sm,n{'^,x,y)\ + \tm,n{u,x,y)\, (4.33) 

we combine (4.26), (4.31) and (4.32) to obtain 

limsup max |/(u, a;, y) -/(a;, y)| < |/(a; - e, y - e) -/(x, y)| 

U—XX) \Tu I Vlog u\/u 

+ lim limsup max (\fm,n{u,x,y)-f{x,y)\ + \sm,n{u,x,y)\ 

«->oo \T^\<Co^/i^/u V 

+ |tm.n('u,a;,y)| + |/m.n(M,x - £, y - c) - f{x - e,y- e)|) 

= l/(a:- e,y- e) - /(x,y)|. 

Since the last term —)> 0 as e | 0, we have completed the proof of the second part of the 
lemma. □ 


Now we are ready to prove the main lemmas in Section 3. 


Proof of Lemma 3.1. Let (/'(a, 6) be the density of (Ari(s„), X 2 (t„))^, i.e., 

1 / la2-2r(|T„|)a6 + 6^ 

i-P(Ki) 

By conditioning and a change of variables, the LHS of (3.1) becomes 


max A'i(s) > u, max A' 2 (t) > u 


(4.34) 






22 


V. Zhou and Y. Xiao 


max -^i(s) > u, max X 2 {t) > u Xi(s„) = u -, 

sGs„+ii-2/“l§ tGt„+«“^/°2T U 

y \ f ^ y \ 

X 2 {tu) = u -- ,u -) u~^dxdy 

uj \ u uJ 

f{u,x,y)^{u,x,y)dxdy, (4.35) 


27rY'T^^r2(|^^u2 


exp 


/R 2 


i + »'(ku|) 

where f{u,x,y) is defined in (4.18) with ^ui‘), Vui') in (4.1), and where 

Hu,x,y) 


:= exp 


1 


“ 2(1 - r{\Tu\))ix + y)- 2r(|Tu|)^) |. 


2(1 - r 2 (|r„|)) V m 2 

Since |^(|'^ii|) — p| —>■ 0 as i/ —>■ cx), it is easy to check that 


ma x 

Vlog u/u 


(j){u, x,y) — e i+p 


0 , as M —>■ oo. 


(4.36) 


Recall Ha{-) in (2.1) and f{x,y) in (4.19). Since ^(•), r]{-) are independent, and 


<«')•'"«") ^(1+") W(ff^) ■■')-1 (rf^) “ ‘ 


,tG 


,tG 






where = means equality of all finite dimensional distributions, we have 
/ f{x,y)e^dxdy 

>P^max^(s) > y(t) > y'^dy 


= / ei+pl 

Jr 


= (1 + pfH^ 




( 1 +p)^ 






i^ + pV 


(4.37) 


By (4.35) and (4.37), to conclude the lemma, it suffices to prove 


lim / max_ 

7r 2 |Tu|<Co-\/IoglI/u 


/(m, X, y)(i){u, X, y) - f{x, y)e i+p 
Firstly, applying Lemma 4.1 together with (4.36), we have 

~ x-\-y 

f{u, X, y)4>{u, X, y) - fix, y)e^ 


dxdy = 0. (4.38) 


max_ 

\ru\<CoV^oEu/u 


—^ 0, as M —>■ oo. 


(4.39) 
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Secondly, as in Ladneva and Piterbarg (2000), we can find an integrable dominating 
function g G L(Sf) such that for u large enough. 


ma x 

\'ru\<Co\/\ogu / u 


f{u, X, y)(j){u, X, y) - /(x, y)e i+p < g{x, y). 


(4.40) 


Therefore, (4.38) follows from the dominated convergence theorem. This finishes the 
proof. □ 


Proof of Lemma 3.2. We first claim that for any compact sets § and T, the identity 
il„(§) + i?„(T) - i?„(S U T) = Pf„(S, T) (4.41) 

holds. Indeed, if we let X = suptgg(x(t) — |t|“) and Y = suptgT(x(^) ~ I^T)) then 

ila(S) + ila(T) - 7J„(S U T) = E(e^) + E(e^) - 
= E(e^l{x<n) +E(e'"l{;f>v}) =iIa(§,T). 

Now let Ti = [0, T]^, T 2 = [mT, (m + 1)T] and T 3 = [nT, (n+ 1)T]. Consider the events 


.4=1 

max -^i(s) > u r 

1 


c = \ 

max X 2 (t) > It > 

1 

LtGtu+u”^''“2Ti J 


B = 

D = 


max -^i(s) > u 

S^Su+U~^^°‘1T2 


max ^ 2 ( 1 ) > u> . 


It is easy to check that the LHS of (3.2) is equal to 


F{Ar\Br\Cr\D) 

= [F{A n c) + F{B n c) - P((A u B) n c)] 

+ [P(^ n D) + F{B no)- P((y4 [jB)n d)] 

- [P(A n (c u D)) + F{B n (c u d)) - P((yi u B) n (c u d))]. ( 4 . 42 ) 


Let R{u) 


(i+p)^ ,,-2 

27r-y/l-p2 


exp 


( i+'dh-l)) 


and ga 


(i+p) 



By Lemma 3.1, we have 


F{A^C) = R(gT)H^, (l + 7i(w)), 

\Qoli,C l J \Qo12,C2 / 

P(B n C) = i?(u)iL„, [J^Vl + q2(w)), 

\Qai,ci / \Qa2,C2 / 

P((A UB)nC) = Riu)H^, ((1 + 73(m)), 

\ QaiiCl J \Qot2,C2 / 


where, for i = 1,2,3, 7 i(u) —S' 0 uniformly w.r.t. r„ satisfying |r„| < Coy/logu/u, as 
oo. These, together with (4.41), imply 


P(A (-1(7)+ F{B n C) - P((A U B) n C) 
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= R{u)H, 


Ti 




0l2 \ i 

Qa2',C2/ \Qoti,ci Qai,ci 


Ti T 2 


( 1 + 0 ( 1 )). 


Similarly, we have 

P(y4 (1D)+ P(B (ID)- P((A LIB) (ID) 


= R{u)H, 


T 3 


H„ 


OL2 \ j 

Qq( 2,C2/ \Qai,Ci Qq(i,Ci 


Ti T 2 


(1 + 0 ( 1 )) 


and 


P(A n (c u D)) + F{B n (c u d)) - p((A u b) n (c u d)) 


= R{u)H, 


Ti UT 3 




012 \ I -^-^Oi 

Qo2,c 2 J \Q(yi,ci Qai,ci 


Ti T 2 


( 1 + 0 ( 1 )). 


By (4.42) - (4.45), we have 

F{AnBncnD) 


= R{u)Hai 


Ti T 2 


^Qai,ci Qoi,ci 

which concludes the lemma. 


H, 


Oi2,C2 


Ti T 3 


Qo2 ,C2 Qa2 ,C2 


(1 + 0 ( 1 )), 


(4.43) 


(4.44) 


(4.45) 


□ 


Proof of Lemma 3.3. Let f{\t\) = yqiyQTiy. Recall r^i defined in (3.9) and Irkij < 
2S(u), when u is large. By Taylor’s expansion, 

/(kkil) = /(O) + i/"(0)|Tkip(l + 7kl(M)), 

where /(O) = /"(O) = as u —>■ 00 , 7ki(rt) converges to zero uniformly 

w.r.t. all (k, 1) G C. Therefore, for any e > 0, we have 

^ ^ (4 4g) 

(k,l)GC (k.l)GC 

when u is large enough. For a > 0, let 

h{u,a):= ^ ( 447 ) 

(k,l)GC 

In order to prove (3.12), it suffices to prove that 

lim di {u)d 2 {u)h{u,a) = ( —) ^ mesN^Ai (^A^). 

«—>oo Va/ 

To this end, we write 

d^{u)d 2 {u)h{u,a) 


(4.48) 
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= ^ y • (udi(w))^(ud2(u))'^. (4.49) 

u , . 

(k,l)GC 


Let 


p{u) 

q{u) 


—^ min e 


max e 

“ (k,l)GC 


(udi {uy{ud 2 iuy, 

{udi{uy {ud2{u)^. 


It follows from (4.49) that 

p{u) < di {u)d 2 {u)h{u,a) < q{u), 

and 


p(«) < ^ <«(»). 


Observe that 


/ I s£uAi^t£uA2 
|t—s| <Cv^log u 


|t—s|<C-\/log 




/ fy^uAi ,x-\-y^uA2 
\x \<C\/\ogu 




e dxdy 


,N 


' \x \<C\/\ogu 


' \ x\ <Cy/TogU 


6 ^ ^ dx I ^{y^uAir\{uA 2 —x)}dy 

Jrn 

6 ^ ^ dx I ^{zGAin{A2 — x/u)}dz 

Jr^ 


mesN(,Aif\ A2) / e dx = (—'j ^ mesAr(Ai O A2), 
jRiv V a / 


(4.50) 

(4.51) 


(4.52) 


as It —>■ 00 , where the convergence holds by the dominated convergence theorem. Indeed, 
/rw ^{z(^ AirM,A 2 -x/u)}dz is bounded by max|£|<i mesN{Ai O (A 2 - e)) uniformly for |x| < 
C Vlog u when u is large enough. 

It follows from (4.50)-(4.50) that, for concluding (4.48), it remains to verify 

D{u) := q{u) — p{u) —>-0, as u ^ 00 . (4.53) 

Define 

P := |(s,t) G Ai X ^2 : |t — s| < d(u) + \/]Vdi(M) + •\/7Vd2(w)|- (4.54) 

By the definition of C in (3.7), we see that V C x Ap^ C V. Since di(u) = 

o(S(u)) and d 2 (u) = o(S(u)) as u —>■ 00 , the set X> is a subset of V := {(s,t) £^ 1 X^ 2 : 
|t — s| < 2 i5(m)} when u is large. 
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Write D{u) in (4.53) as a sum over (k, 1) G C. To estimate the cardinality of C, we 
notice that 


mes2N{'D) = I ^{\t-s\<2S{u)}dsdt 

'seAi,teA2 


' \x\< 2 S(u) J yeAir\(A 2 — x) 


dydx < K S(u 




(4.55) 

(4.56) 


for all u large enough, where K = ^(7V/2) max|j| <imesi\f(Ai n {A 2 — e)). 

Hence, for large u, the number of summands in (4.49) is bounded by 


#{(k,l) I (k,l)GC}< 


mes2N(V) 


< 


K 6{i 


\N 


mes2iv(A^^^ x Ap^) <(u)d^ (u)' 


(4.57) 


Next, by applying the inequality e ^ — e ^ < y — x for y > x > 0 to each summand 
in -D(n), we obtain 


max 

{s,t)£uA^^ 


- min 

{s,t)£uA^^ 


< a 


max 


|t-sp- 




lt- 


= amax (|t — s| + \ti — si|)(|t — s| — |ti — si|). 


(4.58) 


where the last maximum is taken over {s,t,si,ti) S uAp^ x uAp^ x uAp^ x uA[^^ 
Since |t — s| < 2S{u) for all {t,s) G x uAp^ when u is large, the inequality 

||t — s| — |ti — si|| < |t — ti| + |s — si| implies that (4.58) is at most 


(4.59) 


N 


^ay/Nu^5(u)(di{u) + d2{u)) 

when u is large enough. By (4.59) and (4.57), we can verify that 

D{u) < A ay/l^u^5{u){di{u) Y d 2 {u)){udi{u))’^ {ud 2 {u)) 

U d-^ \^Zlj(l2 \dj 

N + 1 2 I _2 _ ^ _ 2_. 

< C'o(logM) 2 [u “1 + M “ 2 ) —>.0, as u —>■ 00 . 

Therefore (4.48) holds. Similarly, we can check that the same statement holds while 
changing the set C to C°. □ 


Proof of Lemma 3.4- Inequality (3.22) holds immediately by Lemma 6.2 in Piterbarg 
(1996). Hence we only consider the case when m 7 ^ 0. Suppose that {X{t),t G M'^} 
is a real valued continuous Gaussian process with E[A(t)] = 0 and covariance function 
r{t) satisfying r{t) = 1 — |t|“ + o(|t|“) for a constant a G (0, 2). Applying Lemma 6.1 in 
Piterbarg (1996), we see that for any S' > 0, 

P( max X(t) > u, max A(t) > u 

VtG«-2/“[0,S]'V tG«-=/“[mS.(m+l)S] 
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max > u I + ] 


max w 

[mS,(m+l)S] 


— Pf max > m) I 

VtG«- 2 /“([ 0 ,S]'^U[mS.(m+l)S]) / 

= (7J„([0, 5]"^) + (m + 1)S]) - H^{[0, Sf U [m5, (m + 1)^])) 

X ^^e-5“'(l + o(l)) 

y/ZTTU 

= iJa([0,5']^, [mg', (m + 1)5']) (1 + o(l)), as m —>■ oo, 

y/ 2 TTU 


(4.60) 


where the last equality holds thanks to (4.41). 

On the other hand, by applying Lemma 6.3 in Piterbarg (1996) and the inequality 
infsg[Q rn+i] |s — t| > l^iol ~ 1 (recall that io is defined in Lemma 3.4), we have 


P( max X(t) > u, max > u 

< exp (-i(|m.J - 1)“5“) 

for all u large enough. It follows from (4.60) and (4.61) that 


iL„([0,5]^,[m5, (m+l)5]) < Coexp 



- 1)“5^ 


(4.61) 


(4.62) 


which implies (3.24) by letting 5 = 

When jmipj = 1, the above upper bound is not sharp. Instead, we derive (3.23) in 
Lemma 3.4 as follows. For concreteness, suppose that io = N and = I- By applying 
Lemmas 6.1 - 6.3 in Piterbarg (1996), we have 


max Xlt) > u, max ^(t) > u 

tGu-2/<»[0,S]«‘ iGii-2/“[mS.(m+l)S] 




max 


tG«“^/“(n(Li''[™j'S’,(m3 + l)S]x [S,S+-\/S]) 


X{t) > 


+ P( max X(t) > u, max ^it) > rt 

VtG^i- = /“[0.S]« tGu-2/<.(nf="d[mjS,(mj + l)S]x[S+VS,2S+\/S]) 

y/27TU v27rii 


< Cn5^-^ 


1 




-e ^ 


(4.63) 


for u and 5 large. Hence, when jmipj = 1, we have 

iL„([0, Sf, [m5, (m + 1)5]) < CqS^-^ 


(4.64) 
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for large S. This implies (3.23) by letting S = 

Notice that 

#{meZ"| max |mi| = fc} = (2fc+1)^ - (2fc-1)^, fc = 1,2,.... (4.65) 

l<i<N 

By (3.23), (3.24) and the fact that /“ dx ^ as T —>• oo, we 

have 


m^O k—1 \miQ\=k 


tXJ 

< Co{3^ - + CoY^[{2k + 1)^ - {2k - 

k^2 

/OO cx a 

dx< CqT^-^ 


for T large enough. This completes the proof of Lemma 3.4. 


□ 


Proof of Lemma 3.5. The proof is similar to that of Lemma 3.3. Indeed, we only need 
to modify (4.49) and (4.55) in the proof of Lemma 3.3. For any y = (yi,...,?/ at) G 
and 1 < J < j < IV, let yi-j = {yi, ■■■,yj). On one hand, with a different scaling, h{u,a) 
in (4.49) has the following asymptotics: 




,x-\-y^uA2 
I cc I < C \/IoglI 


dxdy 


^d^ {u)d^ {u)h{u,a) 

|x|<C 

\ jR" 

IT / ^{Vi^l^Sj,uTj]nluTj-Xj,uRj-Xj]}dyj I dx 

^-a\x\ TT 1 f / 1/ ; .xdzi-M^ dx 

|2l;M6-4i,Mn(A2,M-a;i:M/u)| J 


1 

Tm 


^ -/I a; I < C VTogH 
N 
X 

j^M +1 ’ 


' \x \<C\/\ogu 


c ^ 

TOesM(^i,M n ^2 ,m) / dxi:M IT / Xje~°'^idx 

•_ A/f I 1 


j=M+l ' 


_ 2^ N^M/2^M/2 ^77ieSM(^l,M n A2 ,m), 


(4.66) 


as u —>■ OO. On the other hand, when u is large enough, mes 2 N{P) defined in (4.55) can 
be bounded above by 


mes2N(JA) = 


i^)= II 

J Js^Ax,t^A2 


l{|t-s|<25H} 


dsdt 
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N 

dyi-.M^ n ^V{xj>o}dx 




r r ^ 

= 5{uf^~^ j ( dyi-.n) H ^j'^{^j>o}dz 

d\z\<2 dyi.j^^Ai^MC\(A 2 ,M—Zl:MS{u)) J —M+1 

<KS{u)‘^^-^, (4.67) 


where K = max|e|<i mesuiAi^M n (A 2 ,m - e)) I\z \<2 O^m+i zA{zj>o}dz. 

By (4.66) and (4.67), (3.36) can be obtained through the same argument in the proof 
of Lemma 3.3. We omit the details. □ 


We end this section with the proof of Lemma 4.2. 

Proof of Lemma 4.2. Let /«,ru(-) and /(•) be the density function of X{u, r„) and X, 
respectively. It suffices to prove that for all x G , 


l{y<x} 


f{y) max 




fiy) 


- 1 


dy —>■ 0, as w —>■ oo, 


(4.68) 


where {y < x} = 

First, we will find an upper bound for max,-^ l/u.T„(y)//(y) — 1|- For any e > 0, define 


r(M, Tu) = {'-fij{u,Tu))i,j=l,...n ■= -(S(u,r„)-E) 
e(u, r„) = {ei{u, Tu))i=i,...n ■= ^(m(w, t„) - fi). 

By Assumption (4.4), there exists a constant U > 0 such that for all u > U, 

max|/rj(u, r„) -/ijl < e, max |cry (u, r„) - cr^ | < e, i,j = l,...,n, 

Tu Tu 


which implies |7ij(u,r„)| < 1 and |ei(tt,T„)| < 1 for m > t/. 

Let = {vij)ij=i^,,,^n be the inverse of S. When e is small, the determinant of 
E(u,t„) satisfies 

|E(u,r„)| = |E + er{u,Tu)\ = |E|(1 + etr(E“^r(u, r„)) + 0(e^)), 

where 0(e^)/e^ is uniformly bounded w.r.t. for large u (see, e.g., Magnus and Neudecker 
(2007), p. 169). Hence, when e is small enough, we have 

< 2e|tr(E“^r(w,r„))| < 2e^ \vij\. 

i,3 


\^iu,Tu)\ 

|S| 


(4.69) 
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Since \^ij{u,Tu)\ < 1 , Vi, j = 1 , Vr^ for large u, as e ^ 0 , the inverse of 'Y{u,Tu) 
can be written as 

E(w, T„)-i = E-i - eE-ir(M, r,)S-i + 0 {^), 

where is a matrix whose entries are uniformly bounded and independent of 

for large u (see, e.g., Meyer ( 2000 ), p. 618 ). Hence, 

du,r^iy) ^ {y - h-(u,Tu)Y'Y~^(u,Tu){y - y{u,Tu)) - {y- yfY-^^y - y) 

= -\{y- yf {- eS"^r(M,T„)S-i + 0 {e^)){y - y) 

+ 66 ^( 14 ,T„)(S“^ - eS“^r(it,T„)E“^ + 0 {e^)){y - y) 

- ie^e^(M,T„)(E“^ - eE“^r(M,T„)E“^ + 0(e^))e(u,r„). 

Since | 7 y (it,r„)| and \ei{u,Tu)\ are uniformly bounded by 1 w.r.t. r„ for all u > U,we 
derive that for any y G M^, 


max |d„_T-^(j/)| —>-0, as u —>■ oo. 


(4.70) 


By (4.69) and (4.70), for y G 


max 


f‘U.,T^{v) , 

= max 


f{y) 

T-u 

|E|-V2 


0, as u —>■ 00 . 


If we could further find an integrable function g(jj) on 

fu,T^ (y) 


f{y) max 


/(y) 


- 1 


< y(y), 


(4.71) 


(4.72) 


then (4.68) holds by the dominated convergence theorem. 

Given a constant Co, let Aj := \ maxij \aij\ < Co}, &/ := 

{(^i)r=i € I maxi \bi\ < Co}. Then there exist constants C 2 ,C 3 , such that 

\x'^Ax\ < C2x'^x, \b^x\ < C3 + x'^x, Wx G jWA G Aj,yb G bj. 

Hence, there exists a constant C 4 > 0 such that 

\du,Tuiy)\ <Cie{y - y)'^(y - y) + Cie, (4.73) 

By (4.69) and (4.73), for small e and large u, there exists a constant K such that 


max 


fu,Tuiy) _ , 

/(y) 
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then (4.72) holds and hence we have completed the proof. 


□ 
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